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Abstract 

In a continuing effort to understand divergences which occur when quantum fields are confined by 
bounding surfaces, we investigate local energy densities (and the local energy-momentum tensor) 
in the vicinity of a wall. In this paper, attention is largely confined to a scalar field. If the wall is an 
pg . infinite Dirichlet plane, well known volume and surface divergences are found, which are regulated 

by a temporal point-splitting parameter. If the wall is represented by a linear potential in one 
coordinate z, the divergences are softened. The case of a general wall, described by a potential of 
the form z" for z > is considered. If a > 2, there are no surface divergences, which in any case 
vanish if the conformal stress tensor is employed. Divergences within the wall are also considered. 
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I. INTRODUCTION 

Quantum vacuum energy, or Casimir energy, referring to the quantum energies of fields 
in the presence of material bodies or boundaries, is a mature subject [1]. In the last decade 
there have been tremendous advances in both experiment and theory, so that the so-called 
Lifshitz theory [2] has been confirmed at the 1% level, and theoretically there is now the 
possibility of calculating forces between bodies of practically any shape and constitution. 
Even the effects of finite temperature have now been confirmed [3|. Yet, such effects are 
controversial, and there are many issues that are still unresolved. 

One of the issues that has been controversial almost from the beginnings of the subject is 
that of the Casimir self-energy of an object, as opposed to the energy of interaction between 
two or more rigid objects. For example, in 1968 Boyer calculated the self-energy of a perfectly 
conducting spherical shell of zero thickness and found a surprising repulsive result [4] • This 

result has been confirmed by different techniques by many authors since. Yet within a 

l-'i 
decade, the meaning of this result was profoundly questioned |5j; not only is the meaning of 

self-energy rather obscure, but divergences occur whose omission has resulted in controversy 

up to the present time [6|. Some of these divergences are proportional to the volume, to 

the surface area, and to the corners, so-called Weyl terms, which can be unambiguously 

removed. Curvature divergences are rather more subtle, and the reason Boyer obtained a 

finite result was that the interior and exterior curvature contributions cancel. Situations 







without curvature, such as triangular prisms [7| and tetrahedra [8|], have finite calculable 
self-energies when only the interior contributions are included. 

The above calculations refer to the total energies of the systems. Yet, there is much inter- 
est in local quantum energy densities, or more generally, the vacuum expectation value of the 
stress-energy tensor. There are well-known divergences in these as surfaces are approached 
6|. However, most of the work on this subject has studied perfect boundaries, such as ideal 
conductors or Dirichlet walls. Since there are still issues regarding surface divergences that 
are not well understood, which are particularly relevant when the coupling to gravity is con- 
sidered, in this paper we will consider walls that are modeled by potentials that are "softer" 
than such a perfect wall. (Soft walls were considered earlier, but apparently only for the 



global energy 



9|.) In particular, we will consider massless scalar fields in three dimensions 



in the presence of a potential which depends only on one coordinate, z. We follow Bouas et 



al. [lO| and consider semi-infinite potentials, so that the potential vanishes for z < while 
it is a monomial in z for z > 0. For special cases (Dirichlet, linear, and quadratic poten- 
tials) the energy density may be found explicitly in terms of known functions, but in general 
asymptotics yield the information about the nature of the divergences as the region of the 
potential is approached from the left, as well as the divergences in the region of the potential. 



Unlike many previous investigations of this subject, including some by Milton J6|, [ill], we 
precisely regulate all expressions by inserting a temporal point-splitting. Then precise forms 
of the divergences in terms of the temporal splitting parameter are obtained, which exhibit 
Weyl terms, as well as the nature of the singularity at the boundary z = 0. 



the expectec 

Unlike Ref. 10| we consider a general stress tensor with arbitrary conformal parameter ^; 
that reference considers ^ = 1/4, but we find that for C, = 1/6 the divergences that occur as 
the boundary is approached are removed; in any case, they disappear for a potential higher 
than quadratic. 

II. DIRICHLET WALL 

Consider a massless scalar field in three-dimensional space subject to a Dirichlet wall 

, , i 0, z<0, 

v{z)={ (2.1) 

I oo, z > 0. 
The Green's function, the solution to 

— - V + t; j G{x, x') = Six - x'), (2.2) 



has the form 



G{x, x') = I ^^^^-'^^'-'^e^'-^-^'-'^^ai^^ ^'; -), (2-3) 

where, for z, z' < 0, 

g{z, z', k) = e''''< sinh kz>, (2.4) 

where z^, Zy is the lesser, greater of z and z'. Here 

K^ = kl-uj^, (2.5) 

where we have anticipated making the Euclidean rotation (not just a Wick rotation) 

u -> <, {t - t') -^ IT, (2.6) 



so we may regard /t as positive. 

The energy-momentum tensor for the scalar field is 



(2.7) 



where ^ the conformal parameter, which for the conformal value of ^ = | yields a traceless 
stress tensor, t^x = 0. The connection between the classical causal (Feynman) Green's 
function and the (time-ordered) vacuum expectation values of the fields is 



{<jy{x)<f){x')) = -G{x,x'), 



so the one-loop vacuum energy density of the field is 



(2.8) 
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(tOO) = ±. (509'° + V • V) G( 



lAj > lAj , 



+ i^V^G(x,x) 



(2.9) 



Using the Fourier representation (12.31) we have for the energy density 
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-^^ + ^^ + 11^)* 



[z,z 
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(2.10) 



Here, we have regulated the integral by retaining r as a temporal point-splitting regulator, 
to be set equal to zero at the end of the calculation. We now introduce polar coordinates in 
the C-k volume, so that 

C = KCOs9, |k_L| = fi;sin6' (2-11) 

and then the integral over the regulator term is 

2 . 



d cos 9 e* 



Sm KT. 



Inserting the Green's function (12. 4p into the energy integral (I2.10p yields 



(2.12) 



87r2 



duK 



2kz 



4^-1 



1 d^ 



+ 



1 & 



i2 1 



■ sm KT. 



KT 



(2.13) 



K? Ot'^ J K^ Ot"^ 

The two terms in Eq. (I2.13P consist of a ^-dependent term and a constant. The latter is 
just the bulk energy density arising from the free part of the Green's function, 

. ,. 1 



9o[z,z 



-k\z—z'\ 



2k 



This is evaluated as 
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(2.14) 



(2.15) 



which uses the integral 

fOC 

dx sinx = 1. (2-16) 



oo 



'0 

This resuh is just the well-known volume Weyl term. If \z\ ^ r, we can take r — )■ in the 
remaining term, and we find for the Dirichlet wall 

This is exactly the form found near one of the plates for the two-plate Casimir situation 

n 

lll| . Note that this term vanishes for ^ = 1/6, which suggests that it has no significance, 
disregarding gravity. If we keep the regulator, we can integrate over the whole volume to 
the left of the wall: 

"0 1 /-oo 



dziu — Uq) = — - / dn K^ 



, ^ sm KT cos KT sm KT 
4^ h 2 2 

KT {ktY [ktY 



Sttt^ 



(2.18) 



This uses the evaluations 



dK cos KT = 0, / dKKsmKT = 0, / — sinfi;r = — . (2-19) 

Jo Jo '^ 2 

The result (12.181) is exactly the second Weyl term, expressing the energy per unit area for a 
Dirichlet wall. 

In exactly the same way we can compute all the components of the stress tensor. The 
result, for \z\ ^ r, is exactly as expected: 

{tn = ^;^diag(3, 1,1,1) + ^^diag(-l, 1, 1, 0). (2.20) 

The bulk term has the required traceless, rotationally-invariant form, since it is unaware 
of the wall. The surface-divergent term vanishes for the conformal case, and exhibits no 
force on the wall, so is unobservable. This stress tensor trivially satisfies energy-momentum 
conservation, d^{t'^^) = 0. 

Incidentally, note that if the cutoff were omitted for the bulk term, we would obtain a 
form that is consistent not with rotational symmetry, but with the symmetry for the 2 + 1 
dimensional breakup as seen in the finite part of the stress tensor for the interaction between 
two Dirichlet plates: 



%^-Y^/ dKK'diagil, -1,-1,3). (2.21) 







III. LINEAR WALL 

We next consider the linear wall, 

f 0, z <0. 

v{z)={ . (3.1) 

\ z, z > 0. 

The energy density to the left of the wall is given by Eq. fl2.10p . whereas to the right of the 
wall, the potential must be included, or in general 

u{z) = 1([(9V)' + V0- V0 + t;0V2evV]). (3.2) 

In Eq. (13. 2 p the fields are to be evaluated at coincident points, and again the connection 
with the Green's function is given by Eqs. (12. 8p and (12. 3p . where now the reduced Green's 
function satisfies 

^ +e + v{z)-ujAg{z,z') = 6{z~z'). (3.3) 



dz"^ y 

As we saw before, it is convenient to perform a Euclidean rotation, u —> i(. 

To find the energy density for the region to the left of the wall, z < 0, we solve Eq. (13. 3p 
in the two regions, always assuming z' < 0, in terms of the variable k^ = A;^ + ('^■. 

z<0: g{z, z') = -le""!"-"'! + A(/)e"^ (3.4a) 

2k 

z>0: g{z,z') = B{z')Ai{K'^ + z). (3.4b) 

Here we have chosen the boundary conditions that as 2; — > ±00, the Green's function must 
vanish. The functions A and B are determined by the requirement that the function and 
its derivative must be continuous at 2; = 0. This leads to two equations 

B(z')Ai(K^) = —e^'' + A(z'), (3.5a) 

2k 

-B{z')Ai'{K^) = -—e^'' + A{z'), (3.5b) 

n 2k 

which may be immediately solved: 

^^'^ ~ 2«:" l-A\'{K^)/KAi{Kr ^ ^ ^ 

KZ' 



Thus, in particular, the reduced Green's function in the potential-free region is 



z,z'<0: g{z,z') 



2k 2k l-M{K^)/Kki{K'^y 



(3.7) 



When we insert this into the expression for the energy density fl2.10p we omit the vac- 
uum term in the Green's function, since that has no knowledge of the potential, and was 
completely analyzed in the previous section. We are left with for ^ < (|2;| ^ r) 



u{z) - Uo 



ee 



677^ 



3 2..1 + AiV)A^Ai(^ 

l-Ai'iK^)/KAiiK^y 



(31 



Unlike the integral over real phase shifts 10|] , the integrand is monotonically tending to zero 
as K — )■ oo. The integral is therefore finite for all 2; < 0, and may be very easily evaluated by 
Mathematica. The results are shown in Fig. [H It is seen that the energy density diverges 
as z — )■ 0, not at 2; = 1; in fact, by using the asymptotic expansion of the Airy function. 



l+Ai'{K^)/KAi{K^ 



l-Ai'{K^)/KAi{K^) 8k^' 

it behaves for small negative z like 

1-6^1 
u 



K — )■ CXO, 



(3.9) 



(3.10) 



967r2 z 

The comparison with the exact numerical integration with this leading asymptotic behavior 
is also shown in Figs. [H O 

The solution for the Green's function inside the wall is 

0< z,z' : g{z, z') = 7rAi(fi;^ + z>)Bi{K'^ + z<) - , ,!~ ^]l\'^2 ^^K^^^ + z)Ai{K^ + z') 



(KAi-Ai')(K2 



The energy density within the wall is given by Eq. (13. 2p . or 



(3.11) 



u 



57r^ 



00 /•! 

2 



dKK / dcosO 
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K^ + 2—- + z 



dr^ 



9{z,z) 



d d (9^ 



,JKTCOS 9 



(3.12) 



dz dz' \z'^z g^2 

Because both terms in g involve Airy functions of argument k'^ + z, we can use the differential 
equation for the Airy function to write the above as 



u 



iTC^ 






dKKg{z, z 



sm KT 



(3.13) 




0.000 



FIG. 1: Energy density (divided by 6^ — 1) to the left of a linear potential. The exact result (lower 
curve) is compared with the asymptotic behavior for small z, Eq. (|3.10p . 
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FIG. 2: Relative error of the asymptotic approximation ()3.10p . 



Let us analyze the divergence structure, by considering the first term in g, Eq. (13. lip , 
which would be the term arising if the linear potential existed over all space, because Ai(z) — >■ 
as z -^ oo, while Bi(2;) — > as z — ;■ — oo. In any case, this term corresponds to the bulk 
energy density 

1 , . .X C" . C / " , . , o .^ , n s sm KT 



Uq 



57r^ 



92 92 



dKKnAiiK + z)'Bi(K +z 



(3.14) 
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To see the divergence structure, use the leading asymptotic behavior 

7rAi(«:2 + z)Bi{n^ + z)^ \ S^ , (3.15) 

for large k. Then we write the resulting k integral as 

/ (JK K sin KT = y/z / dy sin ( a/^^ — ly/zr ) , (3.16) 

which for small r is dominated by large y, so that the integral can be approximated by 

The required integrals are, for small /3, 



oo 



1 p^ sin/3y 



rfysin/3y = -, / rfy — p = -/31n/3 + ^(/3), (3.18a) 

1 p Ji y 

/""dyS^^ = -In /3 + constant, H dy""-^^ = ^\n(5 + ff{(5^) (3.18b) 
ii 1/ ii r 2 

and then we see only the r derivative term contributes in Eq. (13.141) . and we obtain the 

expected result [10| 

_ 3 1 z z^ 

27r^ r* 87r^r^ 327r^ 
as the cutoff r — )• 0. (Our point-splittin g pr ocedure would probably not reveal a possible 
5- function contribution suggested in Ref. [lO|.) 



IV. GENERAL z" POTENTIAL 

In general, for an a wall, described by the potential 

f 0, z < 0, 

v{z)=\ (4.1) 

with a > 0, we construct the reduced Green's function in terms of the two independent 
solutions in the region of the potential 



+ k' + z"\ \ ' ' =0, (4.2) 

dz^ J 1 G{z) 

where F{z) is chosen to vanish as z — )■ +oo, and G{z) is an arbitrary independent solution. 
The Wronskian is 

w = F{z)G'{z) - G{z)F'{z), (4.3) 



which is just a constant. 

The Green's function to the left of the wall is 

-'<°. «(-')- 2^-""-"' -^^-"^"ipii!;^- (^■^) 

and to the right of the wall, 

^. ^- > : ,U, .') ^ '-F(..)GM - ^M^M ^jlggj;: (4.5) 

(Adding an arbitrary multiple of F to G, of course, leaves this expression unchanged.) 

For a = 1, F{z) = Ai(K^ + z), G{z) = Bi(/T;^ + z), and w = l/vr, and we recover the result 
in the previous section. For a = 2^F{z) = U{k? /2, \/2z), G{z) = f/(/t^/2, —^/2z), in terms 
of the parabolic cylinder function [12|, |l3| . Alternative notations for this function are 

U{a,x)=D_a-i/2{x). (4.6) 

The value of the parabolic cylinder function, and its derivative, at the origin is 

DM = v^2'^/Vr(l/2-z//2), (4.7a) 

DliO) = -y^e''/^+'/yT{-u/2). (4.7b) 

Therefore, the Wronskian is 

93/2-^2/2 

w = — . (4.8) 

r(fi:V4+l/4)r(fi:74 + 3/4) ^ ' 

The energy density to the left of the wall, 2; < 0, is immediately generalized from Eq. (13. 8p : 
^ ^ " 67r2 Jq F - F'{0) K ^ ' 



For the quadratic wall 

,,-. , 1 2 r{KV4+3/4) 

F(0) + F(0)//t _ l-;: r(.V4+i/4) 

F(0)-F'(0)/fi: 1 I 2 r(KV4+3/4) 

Asymptotically, 

1(^2/4 + 3/4) K / 1 



(4.10) 



1 + 7-1 h «^^> (4-11) 



r(fi:V4 + l/4) 2 V 4fi;V ' 
so we approximate the exact energy density to the left of the wall by 



i-ee rdK i-ee, 

67r2 y^ Sk 487r2 



u{z) -u,^ -^-i -e-^ = --^r(0, -2z\ (4.12) 
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FIG. 3: The lower curve shows the exact energy density for the quadratic wall, for z < 0, using 
Eq. ()4.9p with Eq. (I4.10p . The upper curve is the asymptotic approximation to that energy density, 
given by Eq. (j4.12p . Again the factor 6.^ — 1 is divided out. 

in terms of the incomplete gamma function. The latter is actually a very accurate approxi- 
mation as Fig. |3] shows. 

In the region of the potential, we can calculate the generalization of the "bulk energy" 



Uo 



<,TX^ 






dn kqqIz, z 



smfi;r 



(4.13) 



Q^2 Qri-2 

because the argument leading from Eq. fl3.12p to Eq. flS.lSp holds for an arbitrary potential. 
Here, for the quadratic wall. 



gQ[.z,z 



-UiK^l2, V2z)UUy2, -V2z). 
w 



(4.14) 



The uniform asymptotic approximation for large order for U{k'^ /2, \/2Kt) is given in the 
NIST handbook 12]. The leading approximation is rather immediately found to yield 



1 



9o{z,z) 



+ ... 



(4.15) 



where the subleading term is explained in the following. The leading term differs from 
Eq. (I3.15P simply by changing the potential from z to z"^. This means that we can make the 
same substitution in the integral (I3.17p . and so the bulk energy density (I4.13p is 



3 1 



Uq 



1 z^ 1 

+ 



27r2r4 87r2r2 327r2 



^4 + 2(1 -40-- 



Inr, 



(4.16) 
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where the —(2/3) In r term arising from the subleading term in Eq. f l4.15p results in the 
appearance of the conformal coefficient (4/3) (1 — 6^) for the constant term multiplying Inr. 
This last result may be easily generalized to an arbitrary potential v{z). The bulk Green's 
function at coincident points can be written as 

^ , . 1 



9o[z,z) 



G'{z)/G{z)-F'{z)/F{zy 



(4.17) 



The leading asymptotic behavior of the solutions is given by the WKB approximation 13| . 



F{z) ~ Q-^'\z)e^Y> 
G{z) ~ Q-^/^(2)exp 



dtiq^i^t) 






(4.18a) 
(4.18b) 



where Q{z) = k^ + v{z). Here it was necessary to keep the first subleading correction, as 
given in Ref. 13|. Thus, for large k. 



«''=) ^'<=)-2QV.„)fi +■""(=) 



(4.19) 



G{z) F{z) "^ '~' y ' 8Q^{z) 

This is the immediate generalization of Eqs. (13.151) and (I4.15p . Then, the generalization of 
the Weyl expansion (I3.19P and (14.160 is^ 



Uo 



3 1 



1 V 1 

+ 



27r2r4 87r2r2 327r2 
which uses the evaluation 

duK 



;i-60 






Inr, 



(4.20) 



, -^- — — ^ sm KT ~ -r 

[K^+Vf'^ 6 



=^lnr. 



(4.21) 



which follows from Eq. (I3.18ap . Note that the derivative term vanishes for the conformal 
value of ^. This form, of course, follows from the general heat kernel consideration of this 
problem, and is seen for ^ = 1/4 in Ref. [lO|. 

The behavior of the energy density to the left of the wall is also worked out easily, in 
general. We rescale the equation (14. 2 p for F{z), so that z = x/k, 



dx'^ 



1 + K~2-"a;" F{z) = 



(4.22) 



^ Ref. [lO| proposes that these potential divergences could be subtracted by renormalizing terms in the 
Lagrangian describing the background field v. 



12 



and then solve this equation perturbatively in powers of k ^. F has the form 

F(x)~e-^('l + -i/(x) + ...V «:->oo, (4.23) 

where consistency requires (3 = 2 + a, and / satisfies 

f"{x) - 2f\x) - x" = 0. (4.24) 

This may be immediately solved for /', 



/'(a;) = e'M / dtee-'' + Cy (4.25) 

in terms of a constant C . This will reverse the required decreasing exponential dependence 
seen in Eq. (I4.23P unless 

poo 

C = -l dtee-'^' = -2-^-''T{l + a). (4.26) 

In particular, this determines 







/'(O) = -^±^. (4.27) 

From this we determine the required ratio occurring in Eq. (14. 9 p 



or 



kF{z = 0) F{x = 0) k^+»-' ^ ' 2i+"/t2+"' 



1 + F'(0)/kF(0) r(l + a) 



1 - F'{0)/kF{0) (2fi:)2+" ' ^^-^^^ 

which generalizes Eqs. (13. 9 p and (14. lip . This gives the asymptotic estimate for the energy 
density near the wall on the left: 

1 — 6f 

u{z)-uo —^\z\''-^Til + a)Ti2-a,2\z\), z^O-, (4.30) 

which generalizes Eqs. (I3.10p and (I4.12p . The singularity at 2; = disappears for a > 2; 

"W-'-o-^^T^^-". a>2. (4.31) 



For a <2, 

^^r(l + a) (\zr^Ti2 -a)- - "" 

which as a — !■ 2 from below approaches 

1-6^ 

l8^ 



u{z) -uo = -^^r(l + a) ( |2r-2r(2 -a)- ] , (4.32) 



^(^) - «o = ^^ (7 + ln2|z|) , (4.33) 



an accurate approximation to the general estimate (I4.30p 
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V. CONCLUSIONS 

We have explored in this paper the nature of the divergences that occur in the energy 
density in quantum field theory near walls, for the case of scalar fields. We generalize the 
walls from being perfect Dirichlet boundaries, to potentials of the form z'^ within the region 
of the infinite wall. Besides the usual Weyl volume divergence, which arises from the free 
part of the theory, the energy density exhibits a divergence as the wall is approached if 
the wall is not too soft, a < 2. That divergent term, however, vanishes if the conformal 
stress tensor, characterized by ,^ = 1/6, is used. Correspondingly, there is no observable 
consequence of this surface-divergent term, absent gravity. We also compute the divergences 
that occur within the region of the wall, which depend on the form of the potential. To 
obtain unambiguously observable consequences we would need to consider the interaction 
between two such walls. 

A question arises as to how seriously to take the cutoff. As we noted for the Dirichlet 
wall, if the form for the energy density is taken literally for 2; < r, we obtain a nonvanishing, 
^-independent result for the energy of a single wall, agreeing with the expected area term 
in the Weyl expansion, Eq. fl2.18p . 

How is this analysis generalized for more realistic theories? A similar divergence in the 
energy density occurs near a perfectly conducting boundary if one considers only the electric 
or the magnetic part of the energy, or the TE and TM modes separately [ij]. The results 
there were for parallel plates separated by a distance a, so if we take a — ;■ 00 there we recover 
the energy density for a single conducting wall. The electric and magnetic energy densities 
near such a perfect boundary at 2; = are (the volume divergence is omitted here) 

ue{z) = -um{z) = I ^gM^e-2^^e<^ z > 0, (5.1) 

again keeping the point-splitting regulator. Carrying out the integration, we find 

, , , , f°° dKK^ n^ 2 sin KT 1 T^-12z^ . ^. 

Mz) = -um{z) = 1 -^e __ = _^-,-^-^. (5,2) 

If the regulator is removed, r — )■ 0, 

UE = -UM ^ -^^^^, (5.3) 

the same type of quartic divergence encountered in the nonconformal scalar case, Eq. (I2.17p . 
This result was first observed by Dewitt [15\ more than 35 years ago. Not only does this 

14 



energy cancel when the electric and magnetic terms are combined, but if this energy density 
is integrated over all space to the right of the plate, 

dzuE{z) = —^ I dz = 0, (5.4) 



27r2 7o (r2 + 4^2)3 

we get a vanishing energy! [This result also follows from integrating the integral form of 
Eq. ( I5.2p over z and using Eq. (12. 19^ .] So these surface divergences have but an ephemeral 
existence. (These cancellations do not occur, however, for dielectric interfaces |l4|.) 

If we wish to examine the surface divergences in the complete stress-energy tensor in 
the electromagnetic case, it is better, of course, to break up the description into TE and 
TM modes. For such a local description, we need the rotationally invariant form of the 
electromagnetic Green's dyadic given in Ref. [16i]. Then, it is a straightforward calculation 
using the methods described in this paper to obtain the stress tensor for the TE and TM 
modes in the presence of a perfectly conducting wall at 2; = 0, for |2;| ^ r: 

(^te,tm) = ^diag(3, 1, 1, 1) T ^^diag(2, 1, 1, 0). (5.5) 

Remarkably, both terms have vanishing trace, so the individual modes respect conformal 
symmetry even in the presence of the wall. The 2;-dependent surface term cancels for the 
complete electromagnetic contribution. Neither term would seem to have any observable 
consequence. 
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